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Abstract: By the unitarity cut method, analytic expressions of one- loop coefficients have been given 
in spinor forms. In this paper, we present one-loop coefficients of various bases in Lorentz-invariant 
contraction forms of external momenta. Using these forms, the analytic structure of these coefficients 
becomes manifest. Firstly, coefficients of bases contain only second-type singularities while the ffist-type 
singularities are included inside scalar bases. Secondly, the highest degree of each singularity is correlated 
with the degree of the inner momentum in the numerator. Thirdly, the same singularities will appear 
in different coefficients, thus our explicit results could be used to provide a clear physical picture under 
various limits (such as soft or collinear limits) when combining contributions from all bases. 
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1. Introduction 



In the last ten years, enormous progresses have been made in the computation of scattering ampUtudes 
at both tree-level (see ||l], ^]) and one- loop level (see, for example, the reference |^, ^ and citations in 
the paper). All these progresses come out as a better understanding of the analytic structure of scattering 
amplitudes at various orders. 

At the tree-level, the analytic structure is relatively simple: there is only the single-pole structure. 
However, even with the single-pole structure, since there are many kinematic variables, we are facing 
the multi- variable complex analysis as shown in the old S-matrix program whose central theme is to 
determining scattering amplitudes directly from their analytic structures. This complicated mathematical 
problem is avoided in the BCFW on-shell recursion relations 0, |, |, 0, |l3|, which is an outgrowth 
of Witten's twistor program|14, 15|. The key simplification of BCFW recursion relations is that by a proper 
momentum deformation of two external particles, we have reduced the multi-variable complex analysis to 
the single- variable complex analysis. 

Using the BCFW recursion relations, we can get very compact analytical expressions of tree-level 
amplitudes at the price of introducing the spurious-pole structure. Although each term may contain these 
spurious poles, they will be canceled out after the sum of all terms, since they are not physical. Nevertheless 
they are crucial for simple and compact expressions and have very beautiful geometrical picture as shown 
in |6|,0,|18|,||]. 

At the one- loop level, although the integrand is still rational functions of external momenta (i.e., there 
is only the single-pole structure), its integral will produce the branch cut structure^. The location of 
singularities of one-loop results can be determined by Landau equations Q and these singularities can 
be divided into the first-type and second-type. Based on these analytic structures, a reduction method 
1 20, 21, 22, 23] has been proposed and becomes the standard method for one-loop amplitudes. The 
reduction tells us that a general one-loop scatting amplitude may be expanded in terms of master integrals 
with rational coefficients. This expansion has split the branch cut structure into the master integrals, 
while rational coefficients contain the information of locations of singularities. Since the master integrals 
are relatively well understood ||20| , the one- loop calculation is reduced to the computation of these 
coefficients of the master integrals. Based on the rational structure of one- loop integrands, the very 
powerful OPP method [24, |2^ and Forde's methods 27 1 have been proposed. Based on the branch cut 
structure, a unitarity cut method was initiated by Bern et al [28, 29 1 and was further generalized first in 



36 , 



4 dimensions ||30|, 31 1 and then in generalized D dimensions [32, |3^, 34, 

Although these well known analytic structures of one-loop amplitudes have led us to very powerful 
computation methods such as OPP method and the unitarity cut method, there are still some problems 
regarding the analytic structure to be solved. The first problem is that although from Landau equations we 
can determine the location of singularities, the degree of singularities is not fully discussed. The information 
of degree is very important both for theoretical study and practical calculations, such as the rational part 



^Although the branch cut can be chosen arbitrarily, the starting point (branch-point) has a definite physical meaning. 
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of one-loop amplitudes. The rational part will appear if the expansion master basis is the pure 4D scalar 
box, scalar triangle, scalar bubble or scalar tadpole. As carefully discussed in the rational 

part of one- loop amplitudes contains the double pole like "^j^ ^or ^5(1^, 2+, 3+, 4+, 5+) for the gauge 
theory. Thus to be able to use the recursion relation to calculate, we need to determine the degree of poles 
as well as their residues. Besides the degree of poles, we will meet the structure like which is just 
a phase in the physical region, but a true pole in the general complex plane. The third related analytic 
property we would like to understand is the factorization property when {{K'^ — M^) — )• 0. Unlike the 
tree-level amplitudes with the factorization property A*^^'^ — )• A^£^^A*£^'^, for one-loop amplitudes, we have 
j^i-ioop _^ A^l'^^ A}j^^°°^ + Al^^°°^ A''^^^ + F and there is still no general theory for the structure of the extra 
term T. 

Above discussions are main motivations of our current investigation. To learn more about the analytical 
structure of one- loop amplitudes, it is always very helpful if we can find explicit expressions for one-loop 
amplitudes. With recent developments, especially the unitarity cut methods, now we are able to do so. In 
p^ , analytic expressions for various coefficients of master integrals have been given in the spinor forms. 
However from the spinor forms it is hard to read out analytic properties, thus the first step toward a better 
understanding is to translate the spinor form into the Lorentz-invariant contraction of external momenta. 

After the coefficients are written in manifestly Lorentz-invariant forms, many analytic properties be- 
come obvious. The old analysis made in tells that singularities can be divided into the first- type and 
the second-type. The first-type singularities are fully determined by the dual diagrams and all occur in 
the scalar bases (master integrals) which are well understood. For the second-type singularities, our re- 
sults seem to indicate that they appear only in coefficients of bases. Although it is well known that the 
second-type singularities depend on the dimensionality of space-time, the spins of particles, and the details 
of the their interactions |^], their dependence is not fully discussed. The main reason is that given the 
topology of Feynman diagrams, the denominator of integrand is determined while the numerator can be 
arbitrary. In fact, it is the detail of the numerator that defines the theory under consideration. With our 
general explicit results presented in this paper, we are able to have a better understanding of dependence 
of the analytic structure on numerators. For example, the degree of second-type singularities appearing in 
triangle and bubble coefficients will correlate with the degree of the internal momentum t in the numerator. 
This is one new result coming out from our analysis. 

Having analytic expressions for coefficients will open a door for other analysis. Just like the tree-level 
BCFW recursion relations, we can take pair of momenta to make the deformation, thus reduce the multi- 
variable complex analysis to the single-variable complex analysis. The hope is that with these understand- 
ing, one can find similar recursion relations for one-loop coefficients. Furthermore, the Lorentz-invariant 
forms of one-loop coefficients are also the preparations for two-loop calculations^ using the unitarity cut 
method. 

We must emphasize that although our final goal is to address above problems, the result in this paper 
Recent new techniques for amplitude calculations at two- and higher-loop can be found in B El E| . 
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is just the first step toward this goal and there are still a lot more to be done in future. Thus our results are 
purely theoretical orientated and there is nothing to do with improving the efficiency of current powerful 
methods for one-loop calculations. 

Having above motivations in the mind, in this paper, we show how to transform the spinor form of 
one-loop coefficients into the Lorentz- invariant forms. The evaluation is done within the spinor formalism 



|46|, reviewed in |47]. Using these Lorentz invariant forms, we further discuss the analytic structures of 
coefficients, with some clarifications and interpretations using the S-matrix theory [^]. 

The outline of this paper is as follows. In section 2, we give a brief review of the D-dimensional 
unitarity cut method and the derivation of the one-loop coefficients. At same time, some conventions and 
notations are set up. In section 3, some knowledge about Landau equations and singularities of S-matrix 
programm are reviewed. This section is important to understand our results. In section 4, we transform the 
spinor forms of pentagon and box coefficients into the Lorentz-invariant forms. To enable readers to grasp 
main points of our calculations, we have summarized the result in the first subsection 4.1 and leave the 
details of derivations in later subsections. We do similar organizations in section 5 and 6 for triangle and 
bubble coefficients respectively. Finally in the conclusion section, we summarize main points of the paper 
and have several discussions regarding various possible further applications and clarifications. In Appendix 
A, an typical formula, which is important in the process from the spinor form to the Lorentz-invariant 
form, is given with the proof. 

2. Setup 



In this section, we briefly review the (4 — 2e)-dimensional Unitarity method [B^, 33] and the derivation of 



one-loop coefficients |Q, 35, 36, which are the foundation of our work. In this process, we also set 
up some key conventions and notations for latter calculations. Here we use the QCD convention for the 
square bracket [i j], so that 2ki ■ kj = {i j) [j i]. 

2.1 Unitarity cut method 

The unitarity cut of a one-loop amplitude is its discontinuity across a branch cut in a kinematic region 
selecting a particular momentum channel. By denoting the momentum vector across the cut as K, the 
discontinuity for a double cut can be written as 

C = -^{47rf/' I £^SM{P' " Mf)5W((p - Kf - M|) T{p), (2.1) 

where 

np) = A^I^{p)xA^-f-^^{p). (2.2) 

The T{p) can be calculated by any method, for example Feynman diagrams, ofF-shell recursion relations 
or BCFW on-shell recursion relations M, 01. 
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The "unitarity cut method" combines the unitarity cuts with the famihar PV-reduction method |22|. 
PV-reduction tells us that any one-loop amplitude can be expanded in master integrals li 



A 



l—loop 



(2.3) 



The master integrals in (4 — 2e)-dimensions are tadpoles, bubbles, triangles, boxes and pentagons^. It is 



worth to emphasize that the basis includes the dimensional-shifted basis |51]. More explicitly, if we split 
the (4 — 2e)-dimensional internal momentum p = p/^ + n with p4 the component in the 4D and p the 
component in the (— 2e)-dimension, then the basis include the one like 

,D/2 f d^Pi d-^'t^ if^^y 



IJ^Kp'^Y] can be translated into the scalar basis with dimensional shifting as 



(2.4) 



I^=^-^'[{fi^Y] = -e(l - e)...(r - 1 - e)/^="^-^' "^'lll 



^D=4+2r-2e[ 



(2.5) 



It is, in fact, coefficients of these dimensional shifted bases producing the rational part of one- loop ampli- 
tudes mentioned in the introduction. To see it, let us notice that, for example [52|, 



DiDj 



ITT 

■ 2 

ITT 



d^-^^p 



DiDjDk 
DiD.DkDi 



' 2 , 2 {Pi-Pjf 

rrii + iTLj 



+ 0(e) 



+ Oie) 



ITT 



DiD.Dk 



d p- 



D^DjDkDi 



+-^{Pt+P3+PkT + 0{e) 



(2.6) 



where Di = (p + piY — mf. 

In the unitarity cut method, we derive the coefficient by performing unitarity cuts on both sides of 



Eq.(2.3): 



AA 



l—loop 



(2.7) 



If we can calculate the left-hand side, by comparison, we can read out the wanted coefficients Cj at the 
right-hand side. 



'^For massless external particles, tadpoles do not show up. If we constraint to pure 4D case, pentagons will not show up, 
but rational terms appear. 
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The (4 — 2e)-dimensional Lorentz-invariant phase-space (LIPS) of a double cut is defined by inserting 
two (5-functions representing the cut conditions: 

(2.8) 



(f'^'p 5^+\p^ - Ml)5^+\{p - Kf - M|) 
To simplify LIPS, we can decompose (4 — 2e)-dimensional momentum p as 

p = I+fl- [ (f-'^'p = [ d-'^'n [ d^I, 



(2.9) 



where I belongs to 4-dimensionaI part and /T, (— 2e)-dimensional part. The 4D part momentum i can be 
further decomposed as 



' + zK, 



0; 



d^£= / dz dH6^+\i^) {21 -K) 



(2.10) 



where K is the pure 4D cut momentum and ^ is pure 4D massless momentum, which can be expressed 
with spinor variables as 



f = fP ~ 



P,3; = AA = 1^) 1^] ; J dH5^+\e 
Under this decomposing procedure, Eq.( |2.8| ) becomes 

^4-2e^ = / d^l' ill') 

dfi^ ifi^)-^-' I {£ dl) [I dl] 



dt, {i di\ \ t dt. 



(2.11) 



(4vr)- 

r(-6) 
r(-6) 

where we have used 5-functions to solve parameters t and z as 



2 ^„2^-l-e / ^4^ ^(+)(^2 _ ^2 _ ^^2) ^-(+)((^ _ ^)2 _ ^2 _ ^2) 



1 - 2z)K^ + Ml - M| 



(1 - 2z)K'^ + Ml - M| 



{l\K\l] 



(i^2 + Ml - M|) - v^A[K,Mi,M2] - 4KV^ 
2K2 ^ 



(2.12) 



(2.13) 



with the definition 



/\[K, Ml, M2] = (/C^ - Ml - M'^f - AMIM^ 



-amIm^ 



2M1M2 



2M1M2 
1 



(2.14) 



For convenience, the /i -integral measure can be redefined as 



where the relation between u and /i^ is given by 



f A[K,M,,M2] \ 
V 4i^2 ^ 



U 



It = 



4i^V' 



A[K,Mi,A42]' 



(2.15) 
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Using the new variable u, we can rewrite z, t as 

a — f3\/Y^^ r- K 



z 



where 



a = ^2 ' P = K~2 • (2.17) 



Putting all together, the cut integral Eq.( |2.l| ) is transformed to 



where T{p) should be interpreted as 



T{v) = m , /i^) = ntP,j + zK , = T{\1) , \£] , /i^), (2.19) 



with 

i = tP~ + zK = — — — 

{i\K\i] 



(2.20) 



2.2 Input 

For standard quantum field theory^, T{p) is always a sum of following terms^ 

TU) = . (2.21) 

Y\Ui{^ - K^Y - m'. - ^^^) 

The number of propagators is given by thus to have triangles in the expansion, we need to have k >1. 

To have boxes, k > 2 and pentagons A; > 3. The degree of £ in numerator is given by n + A;, so we require 

n + k >0 only. 

If we define 

n+fc 

i? = ^xji?„ (2.22) 

then 11^=1^(2^ ■ Rj) is just the Oj """^ ^i" component after expanding {2£ ■ i?)""'"'^. So, for simplicity of our 
general discussions, we just need take the following form as the input: 

7-(f) = J2^ • i;)"+^ 



*For non-local theories or some effective theories, the assumption of input in ( |2.2l| ) is not right. 

^Since all external momenta are in pure 4D, the contributions of p in (4 — 2e)-dimension can only have either /i^ -combination 
or £ ■ 7?j -combination. 
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According to the simplified phase-space integration Eq. (|2.18|) , the cut integral can be written as 







mw^^ nil {m\^] 



X / (£ dl) \l dl\ /3Vr^ 'n+2^\ L,^,. - (2-24) 



In the above equation, 



R = |3{^/l^)r + urK, Qi = ^{VT^)qi + aiK (2.25) 



where 



~ R-K R-K 

r = R- -^K, OR = a-^ 

Ki-K K,-K Kf + M?- mf 

<1^ = K^- -j^K, a. = ^-^^ ^ . (2.26) 



For the integrand 



^ = (2.27) 



based on spinor formalism, it can be split into 



where 



= ,.,^:..n+i^; ' > (2-29) 



^ ..s _ V- (A-)"^-(£|Jt;(^,|i;)" ^{t\KR\t} 



The term contributes to pentagon, boxes and triangles, while the Gq term, bubbles. Substituting 
the splitting result into Eq.( p.24 ), and taking the residues of different poles, we can get coefficients of 



various master integrals. Notice that the value of n constrains the basis of master integrals. Terms with 
n < —2 contribute only to boxes and pentagons, and those with n > —1 contribute to triangles in addition. 
If n > 0, contributions to bubbles will kick in. For renormalizable theory, we have n < 2. However, in this 
paper, our discussion adapts to an arbitrary n, such as gravity theory. 
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2.3 Summary of coefficients 

Now we list tlie coefficients of different master integrals. The pentagon and box coefficients are given by 



{Pij,i\^\Pij,2T^'^Y[\=i^tjti,j {Pij,i\Qt\Pij,2] 



+ {Pij,l Pij,2] 



2.31) 



where Pij^i and Pij,2 are two null momenta constructed from Qi and Qj (i < j). More explicitly, if both 
Qi, Qj are massless, then we can set Qi = Pij,i and Qj = Pij,2- If one of them is not massless, for example, 
/ 0, we can construct 



Pij = {Qj + xQi) 

The condition P^j = leads to following two solutions of x 



(2.32) 



Xl,2 



2Qi-Qj±^{2Qi-QjY-AQ}Q] 



(2.33) 



thus we have constructed two null momenta. Formula ( ^.3l| ) makes sense when and only when k > 2. 
Furthermore, if n < —3 (noticing that we need to have n + k >0), there is only pentagon coefficients. 
The triangle coefficient is given by 



C[Qi,K] 



{K 



2\n+l 



2 (VA)-+i(n + l)!(P,,i P,,2)"+' 



+ {Pi,l O Pi,2} 



i^Pi,l-TPi,2 



(2.34) 



where Pj^i and Pi 2 are two null momentum constructed from K and Qi. Formula (|2.34D makes sense when 
and only when k > 1 and n > —1. 

The coefficient of the bubble is the sum of the residues of the poles from the following expression: 



k n 



1 {£\R\ef^^ 



(2.35) 



where poles are given by factors {i\KQi\l) and {i\QtQi\^) ■ Formula (|2.35| ) makes sense when and only 
when fc > and n > 0. 

2.4 Notations 

For convenience, we give notations we will adopt in the paper. First we define the following determinant 



G 



Pi P2 

qi 92 



Pk 
qk 



det {pi ■ qj) 



kxk 



(2.36) 
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If Qi = Pi we write 



G{Pi,P2, ■■■,Pk) = G 



Pi P2 



If Qi =pi for i = l,...,k 



Pi P2 • 

1 in ( |2.36| ), for short we can write it as 

Pi P2 ••• Pk-l Pk 



(2.37) 



{Pk\qk)\ 



Plv,Pfe-l 



G 



Pi P2 ••• Pk^l Qk 



(2.38) 



If the meaning of pi, ...,pk-i is unambiguous, we can even write it as {pk\<lk)- 

Second we define other determinants which are related to the Gram determinant, but depend on the 
masses of propagators. They are 



j^{K,,K,,K;R) 



det 





R-K 
R-Ki 
R-Ki 



K"^ + Ml - M| Kf + Ml - ml K] + M 



K^ 
K-Ki 
K-K] 



m 



2\ 



Ki-K 
K, - Ki 



K,-K 
Kj - Ki 

^1 



with 4 parameters and the structurally similar 





/ K^ 


+ Ml - Ml 


Kf + Ml - 


ml Kj + Ml - 


m| Kl + Ml - 


ml \ 




R-K 


K^ 


Ki-K 




Kt-K 




det 


R-K, 


K,-K 


Kf 


Ki - Kj 


Kt-K, 






R-K, 


Kj-K 


K, - Kj 




Kt-Kj 






\R-Kt 


Kt-K 


Ki-Kt 


Kj-Kt 


Kl 


/ 



with 5 parameters. Another one is 

f K^ + Ml - Ml Kl + Ml 



det 



V 



K-K, 
K-Kj 
K-Kt 
K-K, 



2 
1 

K? 



ml K] + Ml - m] Kl + Ml - mj Kl + Ml 



m 



K, - Kj 
Ki-Kt 
Ki - K.. 



Ki-K 

Kt - Kj 
Kj - Ks 



K,-Kt 



Kt-Kj 

Kl 
Kt-K, 



K,-Ks 
Kj - Ks 
Kt-Ks 
K^ 



(2.39) 



(2.40) 



(2.41) 



which is related to reducing the hexagon to the pentagon. 



3. Singularities 

One main motivation of our calculations it to discuss the analytic structure of coefficients of master inte- 
grals. For this purpose, in this section, we will review some backgrounds coming from the study of S-matrix 
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program ||6|, The main point is that locations of all possible singularities of a Feynman integral can 
be determined, in principle, by the Landau equations. These singularities can be divided into two types: 
the first-type and the second-type. However, as we have mentioned in the introduction, the degree of 
singularities, especially the second-type singularities, has not been discussed in 

3.1 Landau equations 

To start, let us notice that apart from constant multiplicative factors, after Feynman parametrization, the 
general Feynman integral takes the form 



N 



(3.1) 



■0 

with defined by 

N 

ip{p, k,a) = J2 ^ii^li - ^^i)- (3-2) 

i=l 

Here A^, I are, respectively, the numbers of the internal lines and the independent loops of the corresponding 
graph, aj, qi, rrii are, respectively, the Feynman integration parameter, the momentum, and the mass 
associated with the ith line. i'{q) is a polynomial which involves the spins of the participating particles 
and the details of their interactions, n is the dimensionality of Lorentz space. 

The four momentum Qi in any internal line is a linear function of the circulating momenta k and the 
external momenta p. Therefore the quadratic form ip{p, k, a) can be written as 

I I 
ip{p, k,a) = ^ aijkikj + ^ bjkj + c 

i,j=l j=l 

= • Ak-2k^ • Bp + {p^ -Tp- a), (3.3) 

where 

a = aiUi^. (3.4) 

i 

Here, A,B,T are respectively I x 1,1 x {E — 1),{E — 1) x {E — 1) matrices, whose elements are linear 
in a. E denotes the number of external lines of the diagram, k and p are column vectors in the spaces of 
the matrices and their elements are themselves Lorentz four- vectors. 

For most discussions, u{q) = 1 for ( |3.l| ) has been assumed. It is enough for the location of singularities. 
Performing the integration over k in Eq.(|3.lD, we can get 



(3.5) 



DN-{l/2)nl 
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where 



C = det(A), D = -{Bp)'^ ■X{Bp) + {p'^ •rp-a)C, (3.6) 

with X = adj(A)^ and a defined by Eq.( |3.4| ). C is of degree I in the a and D, of degree (/ + 1). Ac- 
cording to the generahzed Hadamard lemma, the necessary conditions for a singularity of / are, using the 
representation Eg . (^^) , 

dD 

Form I : Oj— — = 0, for each i. (3.7) 



If we use the representation (3J.), the Landau equations will be given by 



{ai(qf — mf) = 0, for each propagator i 
(3. 
otiQi = 0, for each loop running by loop-momentum kj 



In both forms ( |3.7D and (3.8), solution with = corresponds to pinch the corresponding propagators. 



so for example, a box diagram will reduce to a triangle diagram. The singularity of a given graph with no 
= (i.e., all propagators are on the mass shell) is called the "leading singularity". 
A connection between these two forms can be found by noticing that an alternative expression for D 
is given by 

D = CD', (3.9) 

where D' is the result of eliminating k from ip by means of the equations 

dtp 

-^j^ = 0, for each j. (3.10) 

In the notation of Eq. ( |3.3D , these equations are 

Ak = Bp. (3.11) 
Together with Eq.(3^) and Eq.( p.lO ), we obtain the Landau equations (|3.8| ). 



OiQi = 0, for each j, (3-12) 

j 

and 

ai{qf — mf) = 0, for each i, (3.13) 
where ^ ■ in Eq.( 3.12| ) denotes summation round the jth closed loop of the diagram. 



X is always well defined even Aet{A) = 0. If det(yl) / 0, we have X = A'^C. 
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3.2 Singularities of the first type 



The Landau equations are usuaUy too comphcated to solve algebraically. So a geometrical method, which 
is the so called dual diagram, have been introduced. The dual diagram is vector diagram for internal and 
external momenta. From dual diagrams we can read out the Landau surface where singularities of the 
first type may locate. For example, for bubble diagram, A[K,Mi,M2] in ( ^.141) is nothing, but exactly 
the Landau surface. From this surface, we can find the location of sing ularities is K'^ = (Mi ± . The 
Landau surface of triangle is given by 



-'tri 



1 -yi2 -yi3 

-y21 1 -2/23 

-ysi -ym 1 



(3.14) 



Pl- 



9 9 

m'^—m'^. 



where yij = yji = '"om^m- ^ with (i, j, k) a permutation of (1, 2, 3). The rrii is the mass of the propagator 
Qi and Pi is the external momentum at the vertex i opposite to the propagator qi. For the box diagram, 
let us denote external momenta clockwise as = Mf, i = 1,2,3, 4 and internal propagators clockwise as 
Qi with mass ( gj-i, q-i and Pi meet at the same vertex), then the Landau surface is given by 



-'box 



1 


-2/12 


-2/13 


-2/14 


-2/21 


1 


-2/23 


-2/24 


-2/31 


-2/32 


1 


-2/34 


-2/41 


-2/42 


-2/43 


1 



(3.15) 



2 2 2 
» 3 



2TOim,- 



where yij = yji 

One important point of the Landau surfaces ( 3.14 ) and ( ^.15[) of the first-type singularities is that they 
depend on masses of inner propagators. 



3.3 Singularities of the second type 

The conventional dual diagrams do not represent all possible solutions of the Landau equations. The extra 
solutions are called the second-type solutions. They correspond to rather special solutions of the Landau 
equations. In Eq.(3.11), if A is non-singular, k will have a unique solution in terms of the p which will 
exactly correspond to the dual diagram construction. Hence second-type solutions will have to correspond 
to A being singular, that is to the condition 



C = det A = 0. 



(3.16) 



Second-type singularities can be divide into two classes, pure second-type and mixed second-type. The 



former, which are given by the Gram determinant equation (2.3( 

G{pi,...,pE-i) = det Pi -pj = 0, z,j = 



l,...,i?-l. 



(3.17) 
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where pi represent any (E — 1) of the E external momenta of the graph. The equation ( ^.17 ) is the 
condition that there be a hnear combination of the vectors pi . . .pE~i equal to zero or, more generally, 
equal to a zero-length vector whose scalar products with pi, . . . ,Pe-i are zero. Detailed analysis reveals 
that second-type singularities stem from super pinches at infinity and correspond to infinite values for some 
of the components of the internal momenta in the Feynman graph. 

Second- type singularities have some properties. First the curve given by ( |3.17 ) is independent of the 
masses of the internal particles. Secondly, the presence of second-type singularities involves the dimension- 
ality of space, the spins of particles, and the details of the their interactions. For example, for pure scalar 
theory, i.e., = 1 in only when E < n {n is the dimension of space-time), second-type singularity 

exists. This result will be changed if ^{q) is nontrivial function. 

In a diagram with several loops, there may be super pinches only for some of the loop momenta 
while the others have ordinary pinches at finite points. These singularities are called mixed second-type 
singularities and their equations will depend upon the internal masses of the lines round the loops with 
finite loops. In this paper, we will focus on one-loop diagrams, so we will not meet the mixed second-type 
singularities. 



4. CoefRcients of pentagon and box 

Starting from this section, we will present the explicit Lorentz-invariant form of external momenta for 
various coefficients of master integrals. Since the transformation from the spinor form to the Lorentz- 
invariant form is a little bit complicated, we will summarize the main results at the beginning of each 
section and leave the derivation in the later part, for which readers can skip safely if they want. 

The pentagon and box coefficients are given by Eq.( p.3l] ). In subsection 4.1, we summarize our results 
and discuss analytic properties of pentagon and box coefficients derived from our calculations. In subsection 
4.2, we discuss how to separate pentagon and box coefficients from the single expression (^) . Then we 
evaluate box coefficients in subsection 4.3. To do so, we need carry out a typical sum, which is done in 
Appendix (see ( [AT| ) and ( [A1T| )). The same sum patten appears also for triangle and bubble coefficients. 



4.1 The summary of main results of current section 

Pentagon: First, the Lorentz-invariant forms of the pentagon coefficients defined by momenta Ki, Kj,Kt,K 
are 



n+k j. G 



Kj Kt K 



/ ^{K,,K„Kt,K;R) V " ^ « \ K, K. Kt / 



where functions G, N can be found in ( ^^Tj ), ( |2^ and ( |2^ . From the expression ( [4.1D following analytic 



properties of pentagon coefficients can be read out: 
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First the factor G{Ki,Kj,Kt,K) is nothing, but the second- type singularity intrinsically related to 
pentagon topology. Furthermore, its degree is {n + k), which is the degree of i in numerator of the 
input (2.21). More explicitly, to have the pentagon in the expansion, we must have /c > 3 in (2.21) 



and n + k > 0. n + k = 0, the singularity G{Ki, Kj, Kt, K) does not appear, but it will be there 
when n + k >\. 

Secondly, there are singularities given by £)i^i<^j<^t,Ku,,K) ^ They come from the trivial reduction 
of the hexagon topology to the pentagon topology and depend on masses of propagators. Their 
dependence of masses is not like that of first-type singularities given in (|2ll ) for the bubble, (|3l4|) 



for the triangle and ( 3.15| ) for the box. In fact, the trivial reduction from the hexagon to the pentagon 



is intrinsically related to the four-dimensional space-time. Thus we guess the appearance of this type 
of singularities is also related to the space-time dimension. 

• Thirdly, C[Qi,Qj,Qt, K] does not contain u at all. In other words, dimensional shifted bases exist 
only for box, triangle, bubble and tadpole topologies. This is tightly related to (4 — 2e)-dimensional 
analysis. In fact, it is well known that if we do reduction in pure 4-dimension, the pentagon will not 
be a basis at all. 

It is worth to point out that since the (^]^) is given by only one term, we will not expect cancelation 
of any factor in denominators^. 



Box: The true box coefficients are given by 



ClQi.Qj.K] 



E E 

zi + ...+Zk+s=n+2 h=0 



G 



n 



Ki Kj K 
Ki Kj R K 



\ 



jq{K„Kj,Kt,K;R) 



\ 



J iG{K„K„K)) 



G{K,,K,,Kt,K) \ G{K, , Kj ,Kt,K) 



with h even 



(4.2) 



where T{h) is defined in Eq.(^) and G, N can be found in ( p7| ), (plO| ). The derivat ion of ( |4.2| ) is given 
in later subsections. 

Now the analytic properties of box coefficients can be read out: 



^It is worth to emphasize that in this paper we will not discuss the cancelation of singularities after summing over contri- 
butions from all bases. For example, it is very possible that the singularity like ]j(i^i'Kj'^t,K^,K) canceled out if we 
sum up contributions from all pentagons. 
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First we notice that although there are first-type and second-type singularities in general, the box 
coefficient contains only second-type singularities. The first-type singularity of the box appears only 
in the box scalar basis^. 

Among all second-type singularities, the one given by G{Ki,Kj,K) is intrinsically related to the box 
topology. Among all terms of ( [4.2| ), there is one and only one term with the highest s = n + 2 and 
all other zt = 0. Thus this term can not be canceled by others and the highest degree of singularity 
G{Ki,Kj,K) is (n + 2). 

It is worth to compare the degree of singularities between the pentagon and the box. The highest 
degree of the pentagon singularity is (n + k) (i.e., the degree of i in the numerator) while that of 
the box singularity is (n + 2). Naively, when we do the reduction, one i in numerator will cancel 
one propagator, thus for box we need to cancel {k — 2) propagators, so the remaining degree of i in 
the numerator is (n + k) — (k — 2) = (n + 2) as we expect. However, the degree of the pentagon 
singularity does not follow the rule. We think the reason is following. The naive observation is based 
on the reduction in 4D. If we do everything in pure 4D, the pentagon will not be a basis as given in 



m 



9 



=Y.^dZ='-''[l] + elt'-'-m (4.3) 
1=1 

where /^"^"^"^[l] is finite when e — )• 0. In fact, it is because we do reduction in (4 — 2e)-dimension, 
the pentagon becomes a necessary basis. Based on the observation, we think the naive observation 
is not applicable to the pentagon topology and each i in the numerator does give a contribution to 
the degree of the singularity. 

The appearance of the second-type singularity G{Ki, Kj, Kt, K) indicates the influence of pentagon 
topologies, which will produce the same box topology when pinching one propagator. Just like the 
previous item, for a given singularity G{Ki, Kj, Kt, K) there is one and only one term inside ( [4.2[) 
with the highest degree (n + 3), thus it can not be canceled by other terms. The highest degree (n + 3) 
of the pole G{Ki, Kj, Kt, K) can be understood by the naive observation, i.e., in the reduction one i 
in the numerator will cancel one propagator. To produce the pentagon topology, we need to get rid 



*It is shown in j40[ that for 4D, the scalar box basis does not contain the second-type singularity in physical sheet. However, 
second-type singularities do appear for scalar triangle and bubble bases in 4D. The general condition is that E < D where E is 
the number of external lines and D, dimensions of space-time. It is worth to emphasize that although second-type singularity 
of box does not show up in the physical sheet, it does show up in other sheet. Thus its understanding is also important for 
the study of analytic properties. 



^With formula (4.3), we can change the choice of basis from 7^^*~^'^[1] to 7^^®~^'^[1]. In this paper, we will choose 
/^^^^^■^[l] as our basis to do the PV-reduction. One of the reason is that with this choice of basis, the pentagon coefficient 
will not depend on u, thus contributions to one-loop rational part will be contained completely in the box, triangle and bubble 
parts. 
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of {k — 3) propagators, thus the degree of i in the numerator becomes {n + k) — (k — 3) = (n + 3). 
Each remaining I will produce one G{Ki,Kj,Kt, K) singularity when pinched to box. 

To see the dimensional shifted basis (which is related to rational part of one-loop amplitudes, see 



(2.4), ( |2.6| )), we need to check the n-dependence part in the numerator. From (4^), all n-dependence 
comes from the factor T{h) and the highest degree of u is [(n + 2)/2]. It is also important to notice 
that each u will be accompanied by a factor G{Ki, Kj, K) (see Eq.( [4.21| ), which will reduce the 
highest degree of the pole G{Ki, Kj, K) for these (rational) terms. 

4.2 The separation of coefficients of pentagon and box 

Since the expression ( |2.3l| ) contains both pentagon and box coefficients, the first step is to separate the 
pentagon coefficient from the box. This separation has been discussed in 1^^. However, since we need 
write them more compactly and systematically and we are dealing with the (4 — 2e)-dimensional massive 



case, which is different from the massless case in [48|, we will give the main steps to write out our results 



and leave some details to be referred to |t48 |. 



Expanding numerator: The first preparation for the separation is to expand r (see ( ^.26| )) in the 
basis Qi, qj,qt as ( remembering r ■ K = 0) 

r = at'^^''-''\t + a^-'^'''-^'\, + , (4.4) 

which is equal to the expansion of R in the basis Ki, Kj, Kt, K because qi ■ K = 0: 

R = af-^^'^-^'^^K, + af-^^'^-^'^^K, + af-^.^^-^-'^)^^ + a^f-^-^^'^^^^ET. (4.5) 



By projecting Eq.( |4.5| ) onto the vectorspace orthogonal to K, we can easily check: 



The Crammer rule gives the solution of Eq. d4. 



G 



Ki ... R ... K 



using the notation ( |2.36| ). The denominator G{Ki, Kj, Kt, K) is nothing, but the second-type singularity 
related to the pentagon determined by momenta K, Ki,Kj,Kt. In other words, it can be considered as the 
"finger print" of the related pentagon. 

Using the expansion Eq. ([4.4|), we have 

{Pi\R\P2] = a^'*-*''^ {Pi\Qt\P2] + /J^^-"-"*'^) {Pi\K\P2] , (4.8) 
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with 



uj=i,j,t 



(4.9) 



where has been given in ( ^.40 ) 



Separating box from pentagon: Having explained how to expand R, now we discuss how to separate 



box from pentagon in (1^. First we give an example like 
expand one R we will get 



{RY 



{RY 



{K){Qt^){Qt^) 
{Rf 



10 



First using K, Ki , Kj , Kt^ to 



+ 



{K) {Qt,) {Qt,) {Qt,){Qt,) {K){Qt,) 

For the first term we expand R using K, Ki , Kj , Kt^ while for the second term we expand the R using 
Ki , Kj , Kt2 , thus we get 



{R) , {R){K) 



{Q. 



+ 



t2} 



+ 



{R) , {R) 



Among these four terms, the last one contributes to the box only. For the first three terms, we expand the 
remaining R and arrive 



ci + 



{QtY) 



+ 



+ 



{KY 



{Qt,) ' {QtY){Qt,) 



+ 



C2 + 



{Q 



t2l 



+ 



{R)_ 

(K) 



The last step is to use Qi, Qj, Qt^, Qt^ to expand K for the fourth term 



{Pi\K\P2] 



{Pl\Qt\P2]{Pl\Qs\P2] aK \{Pl\Qt\P2] {Pl\Qs\P2\ 



a. 



+ 



at 



and we arrive 



ci + 



+ 



+ 



{K) ^ {K) 



+ 



C2 + 



+ 



(K) 



Now we have got the complete splitting. The first, the sixth and the eighth terms contribute to the box 
only. The fourth term contributes to the pentagon {K,Ki,Kj,Kt^) only. The second, the third, the fifth 
and the seventh terms contribute to the pentagon {K,Ki,Kj,Kt2) only. In our splitting, we have carefully 
chosen the way to expand R, so that the contribution to the pentagon {K, Ki, Kj, KtY) appears only once 
while the contribution to the pentagon (K, Ki, Kj , Kt^) appears four times. However, it can be checked 
that the sum of these four terms does produce only one term. 



^"For simplicity we have used such short notation (7?) 
obvious. 



(Pi|-R|P2]- By comparing with (2.31), we hope its meaning is 
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The above splitting can be generalized to arbitrary k and n > —2. First we define 



{Pi\K\p,r+'uUt^^,{Pi\Qt\P2\ 

52[n,/c] = ^P,j-2}. (4.10) 

For the simplest example A; = 3 we will have (for example, t = 3,i = l,j = 2) 

..[,3,.gc.K3,i^..K3,m ,4.n, 

where (4.8) has been used. In the above equation, the first term give the true box coefficient, while the 
second term, the pentagon coefficient. The expression of Cs[n, 3] and F[n, 3] can be obtained by induction 
on n as 

F[n,3] = /3fe*,93;r-)"+3^ ^^^^^^^ ^ „fe*-'?3;r-)^(q„g,,,3;.)n+2-s_ 

For k > A, we use the induction on k and get the complete splitting. 

The contribution to the pentagon part has been given in |48| (with a generalization to the massive 
case) and its evaluation gives (4.1) while the contribution to box part is given by the following sum 

zi+...+Zfc+s=n+2 \t=l,tj^i,j J \\ 1\ \ 2\ / 

where s,zt > and in the sum zi + . . . + Zk + s = n + 2, Zi,Zj should be excluded. This formula is completely 
symmetric on t. 

4.3 Evaluation of box coefficients 

Now we need to evaluate ( [4. 131 ), where the sum ^j^]^^^^ + {Pi — )> ^2}^ appears. This sum is a special 
case of the typical sum defined in ( [A.l| ) and its final expression is given in dATTI) . If we put all Qi ^ K 
andT ^ R in (|A1T|) we will get box coefficients. 

However, there is a technical issue related to the u-dependence (i.e., the ^^-dependence part, which 
indicates the dimensional shifted basis) contained inside the definition of R. To have clear separations of 
fx-dependence, we will expand R smartly. To do so, we construct the vector q^''^^'^ orthogonal to all three 
momenta Ki , Kj , K 

('^o)^"*''' = ^.e.upiqU'^K^ = ^.e.^p^K'^K^K^. (4.14) 
Then we can expand R using Qu Qj,K, (go)^"*'-^ and obtain 

{Pi\R\P2] = pVl^a'^o ''""'"''''^ {Pi\qo\P2] + /Sfe''^^'''*''^) {Pi\K\P2] (4.15) 
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with 



where N and G are defined in (|2!39|) and (^ ). The expression ( [4.15D has the clear u-dependence. 
Putting the expansion ( |4.15| ) back, we have 

\h / ID \ „ \ D ih 



(4.16) 
(4.17) 



{Pi\R\P2r 

{Pi\K\P2Y 



(4.18) 



Summing the above result with the term coming from exchanging Pi and P2 and using the formula ( |A.ll ) 
in the Appendix, we have 

iPVT^)HPi\<io\P2]'' , WVT^Y' {PMPi]" 



{Pl\K\P2] 



+ 



(^2|i^|J^l] 



f 2(2i)''(g2)'>{/32(l_„)[(2g,.g,)2-4<?2g2]+4i^2[„^„^(2,^.g^.)_„2g2_^2g2]}h/2 



[(2gi.g,)2-4q2g2]h 



0, 



, for /i even; 
for h odd. 



(4.19) 



Thus 



{Pl\R\P2Y 
{Pl\K\P2]' 



+ 



{K^G{Ki,Kj,K)Y 



(4.20) 



where^"*^ 

T{h) 



/3\l-u)G(Ki,Kj,K)-K^ 
-K^G{Ki,Kj,R,Ky 



2aiajG 



KK, 
K Ki 



a}G{K,Kj)-a]G{K,Ki) 



h/2 



(4.21) 



Combining Eq.( 4.13| ) and ( ^.2[1| ) the box coefficients are given by 



C[Qi,Qj,K] 



E E 

z-i + ...+Zk+s=n+2 h=0 



G 



n 



K^ K, Kj K 
Ki Kj R K 



jy{Ki,Kj,Kt,K;R) 



t=l,tf^i,j 



G{K,,K,,Kt,K) \ G{K, , K, ,Kt,K) 



\ 



h) {G{K.,K,,K)Y ' "^^^'^^^^ 



(4.22) 



^^Also j3,u,ai,aj have K in denominators. It can be checked that the overall T does not have K in the denominator. 
This is important, because the box coefficient (4.22) will not have as its singularity. 
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where T is defined in Eg. ( 4.21 ). The analysis of the singularity structure of ( 4.22 ) has been given in the 
first subsection. 



5. CoefRcients of triangle 



The triangle coefficient is given in ( |2.34 ) and we recall here that when n > — 1, 



C[Qi,K] 



2 (^/A)"+i (n + l)!(Pi Pa 



d 



'n+l 



{£\RQ,\i) 



n+k 



+ {Pi ^ P2} 



e^Pi-TP2 



(5.1) 



T-S-O 



where Pi,P2 are two null momenta constructed from Qi,K as given in ( 2.32 ). More explicitly Pi and P2 
are given by 

Pi,2 = Qi + xi^2K (5.2) 

with 



3^1,2 = , V A = /3V1 - -uvd, 6 



AqtK 



(5.3) 



To have a good separation of the w-dependence, we can also construct two null momenta pi^2 from qi and 
K 

V6 



pi,2 = qi + yi,2K, yi,2 



"2i^2 



Comparing definitions ( |5.2| ) with ( |5.4| ) we have 

Pi,2 = /3\/l - ^iPl,2, 



(5.4) 
(5.5) 



thus the triangle coefficient ( |5.1D p8[ can be rewritten as 

(K2)"+i 1 1 



C[Qi,K] 



2 (^/^)"+i (n + l)!(pi P2> 



n+l 



-/n+l 



{i\rQi\i) 



n+k 



+ {pi O P2} 



where 



r = r Qi, qt = qt ' 



(5.6) 



(5.7) 



The good property of expression ( |5.6| ) is that only QiS have the u-dependence. Now we will evaluate 
residues based on this expression ( |5.6D . 

The presentation of this section is following. In the first subsection we present the result and analyze 
the singularity structure. For readers who cares only the result, reading this subsection is enough. In 
subsection 5.2, we will evaluate the derivative part in the expression ( |5.6[ ) and finally we give the Lorentz- 
invariant form in subsection 5.3, where the polynomial property of u is a natural by-product. 
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5.1 The summary of main results of current section 

The Lorentz invariant form of external momenta of the triangle coefficient is given by 

{n + k)\ms,s')T2{zt) 



C[Qi,K] 



(K2)2{n+1) g g ^ 



(_2)n+l 



n 



S=Os'=0 {ZI.Z2 zfe}>0 

T,t=l.t^i''t=n+l-s 



2s')\(n + k - s + s')\ 



l+zt 



{1 + 21, 1 + 2:2, .-,1 + 2*} 

E 



T3{zt,ht) 



{hi,h2,---,hk)>0 



G{Ki,K) 



n+l-h/2 



ht=cvcn 



(5.8) 



where h = Ylt=i tj^i Various functions Ti, T3 can be found in ( 5.24| ) and G is the Gram determinant 
defined in ( 121361) and (Jp^ ). 

From ( |5.8| ) we can easily read out the analytic structure of triangle coefficients: 



First the coefficient contains only second-type singularities and the first-type singularity related to 
the triangle topology appears only in the triangle scalar basis (with dimensional shifted basis). 

There are only two kinds of second- type singularities. The first kind of second- type singularities 
is given by G{Ki,K) = 0, which is the second-type singularity intrinsically related to the triangle 
topology specified by momenta Ki,K. The highest degree of the pole G{Ki,K) is (n + l). It fits with 
the naive observation in the reduction, i.e., among (n — k) inner momenta i in the numerator, (k — 1) 
of them have been used to remove {k — 1) propagators, thus it is left with (n + 1) i in numerator 
contributing to the triangle topology. Each i will bring one factor G{Ki,K) in the denominator of 
the coefficient, thus we will have the degree (n + 1). 

For the pole G{Ki, K, Kt), which is the second-type singularity intrinsically related to the box topol- 
ogy specified by momenta Ki,Kt, K, the highest degree is (n+2). Its appearance is very natural since 
these boxes can be reduced to triangled by pinching one propagator. In other words, their influence 
to the triangle is given by the appearance of the factor G{Ki, Kt, K). Moreover, it fits with the naive 
observation in the reduction and is, in fact, the same highest degree found for the box coefficient in 
the previous section. The same highest degree (n + 2) is also necessary for the cancelation of soft or 
collinear singularities between box and triangle contributions. 

It is worth to mention that when k = 1, there is no box coefficient at all. From (5.^ 
the pole G{Ki, K, Kt) will not appear, which is consistent. 



we can see that 



Similarly to the box case, we need to check the u-dependence part in the numerator. From (p.^), all 
u-dependence comes from factors Ti,T2 and its highest degree is [(n + l)/2]. It is also important to 
notice that each u will be accompanied by a factor G{Ki,K) (see Eq.( |5.24| )), which will reduce the 
degree of the pole G{Ki,K) for these (rational) parts. 
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5.2 Evaluation of the derivative part 

The ( |5.6| ) contains the standard sum defined in (A.l), but there is also the differential action. Thus to 
apply the result in Appendix, we need to evaluate the derivative part first. Let us define 

/ = (^|fQ,|^)"+'= , g = i , (5.9) 

then the derivative gives 

^="f .. (5.0) 

where (*)^*'' denote the s-th order derivative of the function (*). The Qi is a linear combination of pi^2 



Qi = l^iPi + l^2P2, Mi,2 = t: — ± (5.11) 

2 lyi 



The evaluation of Z^'*); After some algebraic manipulations, we can easily get 

{pi - Tp2\rQi\pi - TP2) = {pi P2) ao{T - to,i)(t - to,2) (5.12) 

where 



and 



oo = m {P2\r\pi\ , ro,i = ro,2 = ■ (5.13) 



0(r) = (/i2 {P2\r\p2] - fJ'i {Pi\r\pi]f + 4/ii/i2 {P2\r\pi] {pi\r\p2 



4i2r- q^f + 4 (/32(1 - n) - 4 ) " " <ltr') (5.14) 



with the explicit u-dependence. 

To continue, we need the following formula 



{b,b2 . . . 6„)(^) = -r^ — • • • (5.15) 

Zi+Z2 + ...+Zn = k 

After we set 61 = 62 = ... = ao(T — to,i)(t — tq^2)^ to have nonzero result, < zj < 2. Using s' to denote 
the number of 6j having the second-order derivative (so there must be {k — 2s') of bi having the first-order 
derivative), we have 

. . . _ , X 

(6i62...6n)(') = Yl (^,j(^''_2',/J^Kao)^'[ao(-To,i-ro,2)]'=-'^'[aoTo,iTo,2)]"^^^^ (5.16) 
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where we have take r — )• 0. Substituting n ^ n + k,k ^ s into the above result, we obtain 



s/2] 



/(«) = {p, p2>"+' 



n + k\ fn + k — s' 
s - 2s' 



yi 



s-2s' 



-M2 {Pi\r\p2]) 



n+k—s 



(5.17) 



where we have used the short notation (ri|r2) = {ri\r2)qi,K defined in ( |2.38 ) since in this section, qi,K are 
fixed. 



The evaluation of g^"'~^^ Similar to f, g can be written as 

9 = n ^(—^ —) (5-18) 

{pi P2) tJiJ^i " - T - n,2j 

where at, Tt^i andr(^2 are given in ( ^.13| ) with r replaced by qt- Then from Eq.( |5.15| ) we can get 



k 



9 



(n+l-s) 



(n + l-s)! ^ -p-r 1 /I 1 



E n 



zt>0 



:,2 



(n + l-s)! A 1 / 1 1 \ . X 

-^.,S..,-.S.7m[w-^]' ''''' 

zt>0 

Substituting expressions of r^^i and Tt^2 yields 

' zt>0 

The final result of derivation: Putting all together and performing a bit algebraic manipulations, we 
can write Eq.(|5.6|) as 



2 (-2gnn+i z-/ ^ s'\(s-2s')l(n + k-s + s'y. 

zt>0 



xTi{s,s')TM) ^M) , (5.21) 
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where the Lorentz invariant forms of Ti , T2 are 



(5.22) 



The n-dependence is entirely in Ti and T2, thus the polynomial property of u is obvious. 



5.3 The Lorentz-invariant Form 



In i ^lTi ), the sum inside the bracket is the standard one defined in Appendix ( [A.l| ). Thus we can use the 
result ( lA.llD given in the Appendix. First noticing that (pi|(7j|p2] = and (p2|%|pi] = by our construction 
(5.4), {pi\qt\p2] = {pi\Qt\P2] and (pi|r|p2] = (pik|P2] in ( 5.21 ). After some algebraic calculations, Eq.( 5.21 ) 
leads to 



C[Qi,K] 



(i^2)2(n+l) 
(-2)'^+! 



n+1 ls/2] 



EE E 

S=0 s'=0 {z-i,Z2,...,z^)>0 



{n + k)\Ti{s,s')T2{zt) 
s'\{s-2s'y.{n + k-s + s') 



n 



1 



Tz{zt,ht) 



l+zt 



( 

{1 + 21, 1 + 22, -.l + Zfc} 



n+l-h/2 



J 



(5.23) 



where h = 'Ylit=i t^i ^t, G is Gram determinant defined in ( |2.36| ) and (2.37) and 



T2{zt) 



Ti{s,s') = (ni{R)) (2n2{R) 

k [zt/2] 



s-2s' 



t=l,tjti 7t=0 ^ '* ' 



It 



n 



1+zt-ht 



(5.24) 



with 



(5.25) 
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It is worth to mention that it can be checked that Qi has {K'^)^ in the denominator and Q2 has (i^^)^ in 
the denominator. When putting back into ( ^.23 ), the factor from Qi, ^2 will be canceled by the overall 
(j^2-j2(n+i) fac^Qj- other words, will not be a singularity for the triangle coefficient. 

Since our main concern is the highest degrees of poles G{Ki,K) and G{Ki, Kf, K), we will discuss how 
to get this information from ( 5.23| ). For the pole G{Ki,Kt, K) there is one and only one term with highest 
degree in the expression ( F23I) which is given by s = 0, s' = 0, = n — 1 and = for r 7^ i, t. In other 
words, the highest degree of the pole G{Ki, Kt, K) is (n + 2). 

For the pole G{Ki, K), since there are many terms contributing to the highest degree in the expression 
(|5.2g ), it will be a little more complicated and we will use another expression to discuss. Since G{Ki,K) = 
K^qf, we will rewrite ( 5.231 ) using qf. Noting that 



Ti{s,s') 



-Qi(r • qi)y, 



Y[ (i + ^t)("i(9t 



i)Y 



therefore, after removing terms with in the numerator, we have 

{K^f^-^^)a-^\-{r.q,)r+^ -^^^ (n + fc)!2-^-' 

^ (_2)n+i(g2)n+i Z^^Z^ s'\{s-2s')\{n + k-s + s')\ 



{zi,Z2,---,Zfc}>0 \t=l,tj^i 
5^' 



With this explicit form, the highest degree of the pole qf (or the pole G{Ki, K)) is (n + 1). 
6. Coefficients of bubble 

After accomplishing the triangle coefficients, the last thing is to derive the coefficient of the bubble. The 
bubble coefficient is the sum of the residues of the poles from the following expression 

This expression is for A; > 1. For k = 0, the answer is very simple and we write down here^^ 

Cmk=o = 2«+i(n + l) ^^-^^ 



^■^With the definition of (2.2(:) and (2.17), we can see the overaU (K^) dependence, which is the only singularity for the 
case fc = 0. 



- 26 - 



As in previous two sections, we summarize the final result and discuss the analytic property in the 
first subsection. The derivation of the result is given in the next three subsections. In subsection 6.2, 
we present the explicit spinor form after the evaluations of residues. In subsection 6.3, we deal with the 
derivative part and finally in subsection 6.4, we translate the spinor form into the Lorentz-invariant form. 

6.1 The summary of main result of current section 

The Lorentz invariant form of bubble coefficient is given by 

k n n—q s [si/2] 

i=l q=0 s=Max{n-2q-lfl} si=0 s{=0 

( G{K,,R,K) y~''-' ^ 2q^U^-nn-q-s T,{z,)Un,r,) 

'zt>0 



k 



t=l,t^i 



GiKi,Kt,Ky 



( \ 

{1 + 21,1 + ^2,. + rr , , . 

>p T3{zt,ht) 



(6.3) 



where Tq^T^ can be found in ( |6.26| ) and T\^T2^Tj, are defined in Eq.( |5.22| ). In the sum ( |6.3D , we need to 
have (n — g — s + 1 + ri — and h = X^f^^ i^i even number. 

From ( |6.3| ), the analytic property of bubble coefficients can be read out as follows^'^: 

• Like coefficients of the box and triangle, only second-type singularities appear in bubble coefficients 
(remembering the first-type singularity of the bubble is {K'^ — {M\ it M2)^)). There are three second- 
type singularities: G(K^ Ki, Kj) related to the box topology, G{K, Ki) related to the triangle topology 
and related to the bubble topology. 

• The highest degree of the pole is n. is the intrinsic second-type singularity related to the bubble 
topology. Its highest degree fits the naive observation in the reduction: to remove k propagators from 
the denominator we need to reduce the fc's H. in the numerator. 

• The highest degree of the pole G{K,Ki) is (n + 1), which fits with the naive observation in the 
reduction. It is also the same as the pole G{K,Ki) appearing in triangle coefficients. Having the 
same highest degree is reasonable when we consider some soft or collinear limits of full one-loop 
amplitudes after combining all contributions (such as box, triangle and bubble) together. 

• The highest degree of the pole G{K, Ki, Kj) is (n + 1). It is worth to recall that the highest degree 
of the same pole in box and triangle coefficients is (n + 2). This is because to get the infiuence of the 
box to the bubble, one further reduction from the triangle to the bubble is needed. 



^The derivation of the highest degree can be found in subsection 6.4. 
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• To see the dimensional shifted basis (which is related to rational part of one- loop amplitudes (^.41)), 
we need to check the u-dependence part in the numerator. The highest degree of u is [n/2]. 

6.2 Spinor form of the bubble coefficient 

The expression (^) is not yet the spinor form of the coefficient of the bubble. We need to evaluate residues 
of various poles from {£\KQi\i) and {i\QjQi\i). Among these two kinds of poles, the contributions of poles 
from {i\QjQi\i) are zero. To see it, let us start with the typical term 



i,q 



and construct two massless momenta as (see ( |2.33| )) 



where 



yi,2 



-2Qi ■ Qj ± \/aM 



2Q 



(6.4) 

(6.5) 
(6.6) 



Then the residues of the poles 2 



Res(Bi,5)| = ± 



n — q-^-k — l 



1 (p^fmptf 



9+1 



where (+)-sign is for the pole P{ and (— )-sign, for the pole . It is worth to notice that the factor 



¥1 



{i\QjQi\i) appears in both Bi^q and Bj^g up to a minus sign, thus Res(Sj^g) 



,{^,f) = -Res(5j-g)| So 
1,2 ^1,2 



£-i>Pl-rP2 



when we sum up all residues, contributions from {i\QjQi\£) cancel. 

Now we consider poles from {£\KQi\£), which has been carefully discussed in Eq.( |5.2[ ). The residue of 
Pi is given by 

^ . _ (-l)"-'^(j^2)2"+2-nxi - X2)"-'^+^ ci"-^ / {£\R\Pi]'^+' \ 

es[ ,,q)\p, P2)"-«A"-'?+i(g + l)(n-Q)! dr^-i [uUt^^{mtQ^\i) m\Pi]'^' J 

and the residue of P2 

_ (-l)"-g+l(K2)2n+2-g(^^_^^^n-9+l ^n-g / (£|i;Q^|^)"-^+fc-l {£\R\P^]''+^ ' 

{PiP2r-'A--i+\q+l){n-qy. dr-'' [nl^^^^^ {l\Q,Q^\l) {£\K\P,r^\ 

Using the relation (^]^) and the corresponding pi^2 in Eq.(^.4|), similarly to the case of the triangle, the 
above two equations can be simplified as: 

(_l)n-g(^2)n+l 



e^P2-TPi 



Res{Bi^q)\p^ 



' (^|fg.K)"-^+'-' (^|/3(7r^)r + afii^|pi]^+'^ 



^£\K\vx 



i<?+i 



(6.7) 
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and, 



Res(Si,g)|p2 



/3(Vr^) (pi P2r-' Vs''''^\q + l)(n - g)! 



9+1 



(6.8) 



Summing all together we have the coefficient of the bubble in the spinor form 

k n 

^[^] = j;^(Res(i?,,,)|p, +Res(i?,,,)|pJ. 

1=1 q=0 



(6.9) 



6.3 Evaluation of the derivative part 

The spinor form (16.91) is complicated and we need to evaluate the derivation first. To do so, we define 



\n—q+k—l 



nUt^i mQi\i) 



w 



{£\K\i] 



9+1 



(6.10) 



and the derivative is given by 



n—q s 



s=0 si=0 



Si 



(6.11) 



The evaluation of f^^'^^ and g^^ ^'^'>: To get Z*-**^-' and g^^ we can use the result in Subsection 
4.2. If substituting si,n — q + k — 1 for s, n + /c in Eq. ( |5.17] ) we get (/^i,2 is given in ( |5.11| )) 



s[=0 



^i.~m) (piklP2]) 



Si - 2s; 



If substituting s — si for n + 1 — s in Eq. ( ^.20D we get 



(6.12) 



T2izt) 

{piP2)'-\, ^ yr^nti,^.(-^2(Pii5*iP2])^+^' 



(s-si) 
9P2 



0<zt<s—s\ 
-s^(s-si) 



where Ti(si,s;) and T2{zt) are defined by Eq.( |5.22D and Eq.( ^.24D . 



(6.13) 
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The evaluation of if^" For Resi?j^q|p^, w is 

wp, = -^-L_(/j(^/r^)2r . Qi - tP{VT^) {P2\r\pi] + aRV6y+\ (6.14) 
So when s > Max{n — 2*7 — 1, 0} 

{n-g-s) _ (g+1)! (-/3\/r^r-'?-^ 



{2q + l + s-n)\ 

x(/3(^/r^)2r • qi - Tf]{VT^) {P2\r\pi] + aRV6f'^+'+'-'' (palrlpi]""^"^ , (6.15) 
and when s < Max{n — 2q — 1,0}, w^p^ = 0. After setting r — )■ 0, Eq. ( 6.15 ) becomes 

Jn-g-s) ^ (g+1)! (-/3Vr^)"-g-^ 

^■^1 (2g + l + s-n)! (V5)9+i 

X (/3(^/r^)2r • qi + aR^/^)2''+i+^-" (p2|r|pi]"~^-^ , (6.16) 

Similarly for KesBi^g\p^, we have 



Jn-g-s) ^ (g+1)! (-/j^T^)"-'?-^ 

'^^^ {2q + l + s-n)l {-^/6)1+^ 

X (/3(^/r^)2r • qi - aRV6f'^+^+'-'' (pi |r Ips]"""'' • (6.17) 

The final result: Now we want to sum up residues of Pi,P2 of -Bj,<j. Up to a common factor, their 
sum is given by 

{pi\r\p2V^'''^^'''^ 



(/3(^/r^)2r • q, + afi^/^)2«+^+^->^ 



n—q—s 



nt=i,^.((pii9*iP2])i+^* 



+(-l)"+^(/3(Vr^)2r • q, - a^V^)2«+i+->';^-'^-,^^^^iL^iL^^^. (6.18) 



nt=i,i^i((^'2igtipi]) 

Using the binomial expansion 

(af n \ s^n-q—s 2q+\+s—nn-q—s 

^^4^- It) = ^ E ^--(^i) (6-19) 
r-i=0 r2=0 



where 



a,..(Pi) = + ^ + ' " ") ("^ (/3(^/r^)2r . g.)'^^(«i?V^)^''+^+^-'^-'^^ 



and similar expression for 

a„.,(P2) = (-i)"-'^-^-^+^^+'^^a,,,,(Pi), (6.21) 
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Eq.( 6.18 ) becomes 



2q+l+s—n n—q—s 

C'ri,r2(pi) 

ri=0 r2=0 



(pik|p2] 



s+k—l—s\ 



J{t=l,t^i{'<Pl\<lt\P2]) 



1+zt 



+ 



(-1) 



n—q—s—l+ri+r2 



{P2\r\pi\ 



\{t=i,t^i{{p2\qt\pi]) 



1 + Zt 



(6.22) 



Since the bubble coefficient is the polynomial of n, from the factor {/3\/T' 



^•^n~q-s-l+n+r2 ^^IG SUm, 

only the terms with (n — q — s — 1 + ri + being even numbers are left. In other words, the above 
expression can be written as 



2q+l+s— n n—q—s 
ri=0 r2=0 



{Pi\r\p2] 



i+k—l—si 



+ 



{P2\r\pi] 



s+k—l—si 



for which we can apply the general expression in Appendix (^) and ( [AllD . 
6.4 The Lorentz invariant form 

Putting all together, the coefficient of the bubble is given by 

k n n—q s [si/2] 

i=l q=0 s=Max{n-2q-l,0} si=0 s[=0 



(6.23) 



G{K,,R,K) 



2 n 



n-q-s 2q+l+s-nn-q-i 

Y Y 

2i>0 



E 

t>i 
(. 



T2izt)T4^{ri,r2) 



G{Ki,Kt,K) 



l+zt 



{1 + 21, 1 + ^2, + 

E 

{h^,h2,...,hk}>0 
\ E*=i,t^i'it=even 



{-AG{Ki, iv:))("-9-^+l+n-r2)/2 



T3{zt,ht 



G{Ki,K) 



n—q—h/2 



(6.24) 



where 



To{s,si,s'i) 
74(ri,r2) 



2"g!(n-g + A;- 1)! 



(2g + 1 + s - n)\s[\{si - 2s[y.{n -q + k-l-si+ s[)\{n - q - s)\ 
la + 1 + s — n\ / n — a — s\ I 2 „ ( K Ki \ \ 



(6.25) 



f 2q + 1 + s — n\ /n — q — s 
ri J\r2 



\K R 



- n))t(""'?-^-i+''i+^2)a''+'+"~"-'^i(-ai) 



1 

n—q—s—r2 



(6.26) 



and Ti,T2,T3 are defined in Eq.(5.22). In the sum ( |6.24 ), we need to have {n — q — s + l + ri— r2) and 
h = X]t=i t^i to be even number. 

From ( |6.24D we can see various second-type singularities. Now we discuss their highest degree. For the 
pole G{Ki,Kt, K), there is only one term to contribute. By setting s = n — g, g = 0, si = 0, ri = 0, 1, r2 = 
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we find the highest degree is (n + 1). It is different from the highest degree (n + 2) of the same pole in 
coefficients of the box and triangle. 

For the highest degree of the pole K^, which is the intrinsic second- type singularity related to bubble 
topology, we can find 



cm ^ jj^ (6.27) 



by noticing 



(^2)2.1 ^ (K2)2^t 

1 



(^J^2'^r-i_+r2 + (n—q-s-l+ri+r2)+2{2qJrl+s-n—ri+n—q—s—r2) 

For the pole G{Ki, K) = K'^qf, its highest degree shows up in many terms, so we need to rewrite the 



result to see clearly. Using 6 = —AqfK'^ and removing all qf terms in the numerator we find 

<5" 



C[K] ^ ^ '\J, ^ «r E 2^'~''^n{s,s„s[){-r-q,r 

s=Max{n-l,0} si=0 s[=0 




From this expression, we can find the highest degree of pole G{Ki,K) = K qf is (n + 1). 
7. Conclusion 

In this paper, to prepare the study of analytic properties of one-loop amplitudes, we have rewritten the 



spinor forms of one-loop coefficients given in |34] to manifestly Lorentz-invariant contraction forms of 
external momenta. Although the rewriting is a little bit complicated and some skills within the spinor 
formalism are needed, the final results of various coefficients are manageable and have been summarized 
in the first subsection of section 4,5,6. 

The main results of our calculations are following. Firstly we have found that although there are two 
types of singularities by the general S-matrix analysis, coefficients of each basis contain only second-type 
singularities. Secondly, the degree of each second-type singularities is tightly related to the degree of 
the inner momentum in numerators. In other words, when we study the analytic property, not only the 
structure of the denominator, but also the structure of the numerator, play an important role. For the 
renormalizable theory, there is an up-bound for the degree of i in the numerator, thus the possible highest 
degree for each second-type singularity is known. For the non-renormalizable theory, the information of 
the degree of i in numerator will also tell us how bad the contributions from these singularities could be. 



- 32 - 



Thirdly, for a given basis, its coefficient contains not only the second-type singularity related to its 
topology, but also those related to its mother topology. Not only that, the degree also matches up. For 
example, the triangle coefficient contains the second-type singularity of the box topology with the same 
highest degree (n + 2) as the box coefficient. This matching has the physically meaningful cancelation 
in various singular limits. To have a clear picture, we have given a table (see Table ||) where for each 
coefficient, the involved singularities and their highest degrees have been given. 



Table 1: The table of singularities and their highest degrees for each coefficient 





j^{Ki,Kj,Kt,Ku,,K) 




G{Ki,Kj,K) 


G{K,,K) 




Pentagon 


I 


{n + k) 








Box 




(n + 3) 


(n + 2) 






Triangle 






(n + 2) 


(n + 1) 




Bubble 






(n + 1) 


(n + 1) 


n 



From the table, there is one point we want to mention. The pole G{Ki,Kj,Kf,K) related to the 
pentagon topology shows up only in the coefficients of the pentagon and box, while the pole G{Ki, Kj,K) 
related to the box topology shows up also in the coefficients of the box, triangle and bubble. We believe 
the reason is that in the (4 — 2e)-dimension, the pentagon can be expressed as the linear combination of 
boxes plus terms in the higher order of e. Thus the influence of pentagon can not be propagated to lower 
topologies. It is consistent with the well known fact that second-type singularity depends on the dimension 
of space-time and structures of interactions, such as the spins, derivative interactions etc. 

We want to emphasize following points for our work. Firstly our paper is to set up a frame for the 
theoretical study of analytic property of one-loop amplitudes. Thus although it will be possible to use these 
explicit Lorentz-invariant forms of coefficients in real one-loop calculations numerically or analytically, we 
will not do it here and leave it as a future project. To use our result in real calculations, we need to discuss 
the efficiency or stability of calculations as carefully discussed, for example, in paper |53]. 

One possible application of our results is following. Since our results are complete, i.e., there are 
(/i^)"-terms corresponding to rational part as mentioned in section 2, we could use our result to calculate 



the rational terms and compare with results from the recursion relation given in [37, 38, 39 1. We believe 
this calculation will help to us to clarify some points in the recursion relation. 

Secondly, our current focus is coefficients of various basis, not the whole structure of one-loop ampli- 
tudes. For the latter, one need to address cancelations of various singularities under various limits (such 
as soft and collinear limits) by combining contributions from various basis. The cancelation of these sin- 
gularities in a complete loop amplitudes is very intricate. It reflects many important information of the 
theory under consideration, such as how good (or bad) the divergence will be and if the theory has some 
unexpected hidden symmetries. It is definitely an important issue and relates to the application in real 
processes as mentioned in previous paragraph. Although our results in this paper provide a starting point 
for these discussions, its explicit demonstration will be very complicated and deserves to be an independent 
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work. 

Thirdly, for coefficients of boxes, triangles and bubbles, there are many terms, thus there are many 
different ways to write down the sum. However, we want to emphasize although there are many different 
ways to group numerators, the denominators are the same. In other words, the appearance of various 
second-type singularities is common for all different expressions, especially the highest degrees of second- 
type singularities are the same. The choice we made here is because we believe this choice gives the best 
presentation of singularity structure. 

Fourthly, as we have mentioned again and again, the coefficients of various bases contain only second- 
type singularities as classified in [^, One exception is the singularities for the pentagon as given in 
jjKi,K,,Kt,K^,K ^ Eq.(|]l])). Although it contains the mass, we do not think it belongs to the first- type 
singularity. This exception is related to the special position of pentagon as a basis in the PV-reduction as 
we have discussed many times in the paper. What it means or if it is a really a singularity deserves further 
study. 

Fifthly, although in this paper, we have identified all singularities. It is still have a lot of work to do 
to understand their properties. For example, we need to know if they are true singularities for one-loop 
amplitudes when we sum all together. If they are, where are their locations: on the physical sheet or 
unphysical sheet. Here we want to distinguish one thing. One of our motivation of current calculations is 
to find a recursive way to calculate coefficients of various basis. Thus all singularities we found in the paper 
do contribute coefficients no matter whether they are physical or not or which sheet they locate at. The 
situation is different from BCFW on-shell recursion relation for tree-level amplitudes where we calculate 
the complete (global) tree amplitudes thus spurious poles do not give contributions. 

Finally, another possible application (which is one of our main motivations) of our results should 
be mentioned. With Lorentz-invariant forms of coefficients, we can study their factorization property 
under the various deformation (such as the BCFW-deformation or Risager's deformation [^]). Based on 
information under deformation, we can try if it is possible to establish some sort of recursion relation. 
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A. Sum in spinor form to Lorentz form 

In this appendix, we will present a formula which is very important to transform the sum in the spinor 
form to the Lorentz-invariant form. The typical sum we meet again and again is the following 

_ {P,\T\P,f , {P2\T\P,f 

= Kf Kf (A-lj 

Uli{Pi\Qt\P2] Uli{P2\Qt\Pi] 
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and 

where Pi and P2 are two null momenta constructed from Qi and Qj. Furthermore we suppose i and j 
are not in the set {1, . . . , N} of ( [A.!]) . Our derivation of the Lorentz-invariant form will use the inductive 
method. 

A.l Recursion relation 

For a given pair (n, m) simple calculations from ( |A.2D give 

vin^y \n^ y ^ (Pi|r|P2] {P2\T\P,] {P2\T\P,f-^ 



{Pi\Qn\P2] (P2ig™iPi] nr=w„ {P2\Qt\Pi] 

^ {P2\T\Pi]{Pi\T\P2]{Pl\T\P2f-^ 
{P2\Qn\Pl] {Pl\Qm\P2] Yit'l^n {Pl\Qt\P2] 

ytn^y rn 1 y ^ (PilTlP^] (P2ITIP1] (P2|T|Pir-^ 

{Pl\Qm\P2] {P2\Qn\Pl] Ul^t^m (^slQil^l] 
^ (P2|r|Pi](Pi|T|P2](Pi|r|P2]^-^ 
(P2|Q^™|Pl] (Pl|Q„|P2] UZu^m {Pl\Qt\P2] ■ 

The sum of Eq. ( A.3 ) and Eq. (1A.4D yields 

^l{Qn)^N-l[Qn] +^l{Qm)^N-l[Qm] = ^Hat + ( ,p ill ip^j/D^iL^ I p 1 ^ '^^ I 2]( 2I I l] 



(A.3) 



(A.4) 



{Pl\Qn\P2] {P2\Qm\Pl] {Pl\Qm\P2] {P2\Qn\Pl] 
X ^N~2[Qn,Qm]- 

Using the spinor formulism (remembering Pi,P2 are two null momenta constructed from Qi,Qj), after 
some trivial manipulations we can get 

{Pl\Qn\P2] {P2\Qm\Pl] + (Pl|Qm|P2] (P2|Qn|Pl] = -^iQn\Qn.) (A.5) 

4 

a 

where we have defined 

^ Qi Qi ' Qj Qi ' Qn ^ 

{Qn\Qm) = det 



(P1ITIP2] (P2ITIP1] = -7^(T|r) (A.6) 



Qi ' Qj Qj Qj ' Qn 

\ Qi ' Qm Qj ' Qm Qn ' Qm / 



(A.7) 



So we get the following relation 



^l{Qn)^N~l[Qn] +'^l{Qm)^N-l[Qm] = 2I!Af + 2 ^ in^Vn^lri \ ^^-2[Q"oQm] (^-8) 
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Summing over all pairs (n, m) of ( [A.q ) we get 



N 



(AT _ 1) ^ Si(Qt)Sjv-i[Qt] = NiN - + 2 ^ 



t=i 



l<t<k<N 



{T\T){Qk\Qt) 
{Qk\Qk){Qt\Qt 



-^N-2[Qk,Qt] (A.9) 



or 



1 



T E 

l<t<fc<Af 



{T\T)iQk\Qt) 

iQk\Qk)iQt\Qt] 



[Qki Qt 



(A.IO) 



A. 2 Proof by inductive method 

With some calculations, we find the explicit expression of S^r to be 



1 



N 



^ ,2^TT(r|Qj) 



[N/2] 

E 

m— 1 



(_l)™2^-'"m!(iV - 2my.AN,,n{T\Ty 
iV! 



N 



(A.ll) 



where the notation means to take the maximum integer equal to or less than N/2, and 



^n-l,m + ^n~2,m-l; 2m < n 

2, 2m = n 

0, 2m > n 



(A.12) 



The second sum at the second line of (A.ll) is over all different choices of m pairs in the set {1,2,..., N} 
and each pair contributes a factor {Qp^ Qp2)- After m pairs having been chosen, each remaining element 
will contribute a factor [Qq T). First few examples A = 1,2,3 can be calculated directly as 



(P1IQ1IP2] {P2\Ql\Pl 

. (r|Qi)(r|Q2) 



{Qi\Qi) 



(T|r)(Q2lQi 



{Qi\Qi){Q2\Q2) (Qi|Qi)(Q2|Q2)' 



23(r|Qi)(r|g2)(T|Q3) - 2(r|r)(g2|Qi)(T|Q3 



(Qi|Qi)(Q2|Q2)(Q3|Q3) 
-2(r|r)(Q3lQi)(T|Q2) - 2(r|r)(Q3|Q2)(T|Qi)), 



which are the same as given by (A.ll). 

We will prove the formula ( A.ll ) by showing that it satisfies the relation Eq. ( A.10| ) by inductive 
method. We check this term by term. For the first term of Eq.( [A.ll ), it satisfies the relation Eq. ( |A.10 ) 
obviously since only the first term of Eq. ( |A.10 ) contributes. For the second part of the formula (A.ll) 
with given m pairs in the set {1,2, ...,A}, both terms of Eq. ( A.lOj ) will contribute. To simplify our 
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discussion, we use the set M = {m pairs} and the set Q = {N} — Ai. The contribution from the first 
term of Eq. ( |A.10 ) is given by 



The sum in the first hne of Ti comes from choosing which q £ Q belongs to the Si part. The contribution 
from the first term of Eq. ( [A.IOD is given by 

i2A^-™-i(m - 1)!(A^ - 2m)lAN-2,m-i{T\Ty 



l[{Qp,\Qp,){N-2m)ll{Q,\T). 



N{N-l)UtiiQk\Qk){N-2y. '.^^'"'""^'psM-p 
-l)'"-i2^-'"-i(m - l)!(iV - 2myAN-2,m-i(T\Ty 



N{N-l)UtiiQk\Qk){N-2) 



The sum in the first hne of T2 comes from choosing which pair p £ M does not belong to the S7V-2 part. 
Summing Ti and T2, with a little algebra, we can see that it reproduces the corresponding terms of formula 



(A^). 



A special case of the above proof is that when N = 2m, only the second term of Eq. ( [A. 10 ) contributes. 
It is easy to see that we do have ^2m,m = ^2m-2,m-i = 2 as given by ( A.12| ). 
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